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ABSTRACT 

In this paper we discuss the contribution of different astrophysical sources to the ionization 
of neutral hydrogen at different redshifts. We critically revise the arguments in favour/against 
a substantial contribution of Active Galactic Nuclei (AGNs) and/or Lyman Break Galaxies 
(LBGs) to the reionization of the Universe at z > 5. We consider extrapolations of the high-z 
QSO and LEG luminosity functions and their redshift evolution as well as indirect constraints 
on the space density of lower luminosity Active Galactic Nuclei based on the galaxy stellar 
mass function. Since the hypothesis of a reionization due to LBGs alone requires a significant 
contribution of faint dwarf galaxies and a LyC photon escape fraction (/esc) of the order of 
~ 20%, in tension with present observational constraints, we examine under which hypoth- 
esis AGNs and LBGs may provide a combined relevant contribution to the reionization. We 
show that a relatively steep faint-end of the AGN luminosity function, consistent with present 
constraints, provides a relevant (although sub-dominant) contribution, thus allowing us to re- 
cover the required ionizing photon rates with /esc ~ 5% up to z ~ 7. At higher redshifts, we 
test the case for a luminosity-dependent /esc scenario and we conclude that, if the observed 
LBGs are indeed characterized by very low /esc, values of the order of fcsc ^ 70% are needed 
for objects below our detection threshold, for this galaxy population to provide a substantial 
contribution to reionization. Clearly, the study of the properties of faint sources (both AGNs 
and LBGs) is crucial. 

Key words: cosmology: observation - early Universe - quasars: general - galaxies: active - 
galaxies: evolution 



1 INTRODUCTION 

, Cosmic reionization is a major focus in present Cosmology, since 
it represents a crucial cosmic epoch for the formation of the first 
structures and the production of the photons responsible for the 
end the Dark post-recombination Ages. Moreover, these energetic 
photons affect, in a critical interplay, the species available for gas 
cooling (and consequently the star formation) and the collapse of 
(small) dark matter halos. 

An important constraint on the epoch when this phase transi- 
tion occurs is set by the measurement of the Thomson optical depth 
of the intergalactic medium (IGM) via the large scale polariza- 
tion o f the Cosmic Microwave Background (CMB. IXomatsu et alj 
[201 l| ), which provides - with the rough assumption of instanta- 
neous reionization - Zreion = 10.6 ± 1.2. Additional evidence 
comes from the Gunn-Peterson test applied at the Lyman-Q for- 
est, chara cterized by a lo w neutral hy drogen fraction at red shifts 
below 6 jpan et aljHooe! , but see also iMcGreer et al]|201 ll) . Re- 
cent evidence for a damping wing around the systemic redshift of 



a z=7.085 quasar jBolton et alj|201 11) and for a sudden decrease in 
the fraction of Lyma n-Q emitters among z ~ 7 Lyman-break galax- 
ies (LBG, see e.g., IPentericci et alj I2OI H; ISchenkeretallbod : 
lOno et aill2012l) also suggests a rapid increase of the neutral frac- 
tion of hydrogen in the Universe at these epochs. All these evi- 
dences broadly constrain the epoch of hydrogen reionization in the 
redshift range 6 < z < 12. 

At the same time, many studies have addressed the nature of 
the ionizing sources. Pop III stars (e.g. [ cia rdi et al. 2000), star- 
forming galaxies (e.g. Robertson et al. 2010l) . and Active Galactic 
Nuclei (AGNs, e.g.|H aiman & Loeblll99i) have found to be prime 
candidates, with more exotic possibilities exp lored in the form o f 
primordial black-holes and mini-quasars (e.g. iMadau et al.l 2004 ). 
and decaying particle s (dark matter and neutrinos, e.g. IScott et al.l 
ll99ll : |PierDaolill2004l) . It is commonplace that by 2: ~ 6 the ion- 
izing radiation em itted by quasars alone is i nsufficient to reion - 
ize the IGM (e.g. ISchirber & Bullockl I2OO3I : ICowie et al.ll2009t) : 
on the other hand high-redshift galaxies are in principle able to 
produce the bulk of the cosmic emissivity ionizing the IGM (see 
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iHaardt & Madaiil ( 1201 2h for a recent analysis), but only if the frac- 
tion of the 1-4 Ryd photon s escaping the galaxi es is significant 
(~ 20% at 2 ~ 7, see also iBouwens et alj 201 If) and if the con- 
tribution of very faint, undetected objects is taken into account, 
assuming their space density is satisfactory described by the ex- 
trapolation of the faint-end of observed luminosity function (LF). 
Nonetheless, the direct detection of Lyman continuum (LyC) pho- 
tons fro m low redshift {z < 1.5) galaxies h as been so far unsuc- 
cessful l lCowie et al.ll2009l : ISiana et alj2010l) . down to fesc < 1%, 
and the measurements at the highest redshifts for which the direct 
LyC measurement is still allowed by the rapidly increasing IGM 
opacity (z — 3 — 4) have been shown to be prone to a si gnificant de- 
gree of conta mination by lower-redshift interlopers (IVanzella et alJ 
I2010al , l2012h . 

The results o f numerical hydrodynamical simulations (see e.g. 
ICiardi et"al]|201 ih also favour sources with a soft spectral energy 
distribution (indicative of Population II stars) and high escape frac- 
tions as the main contributors to Hydrogen reionization. On the 
other hand, numerical simulations have not provided yet a con- 
cordant answer to the question of fesc evolution as a function of 
galaxy properties, luminosities and redshift: some groups reported 
evid ence for a decrease of /esc with decreasing halo mass at 2 > 3 
(e.g. lOnedin et al.ll2008h . while competing groups have found the 
exactly op posite result of an increase of /esc with decreasing halo 
mass (e.g. lYaiimaetalhoul with a large scatter in the individual 
/esc determinations). 

The magnitude of the faintest dwarf galaxies represents an- 
other critical point in the analysis. The expected absolute magni- 
tude at 145 nm (A/uv) of a galaxy hosted in a hal o with virial 
tempe rature ~ 2 x 10* A" is A/uv '-^ — 10 (see e.g. iTrenti et al.l 
I2OIOI) . However, theoretical models also suggest that star formation 
might be strongly s uppressed in early Afuy ^ — 13 dwarfs due t o 
metallicity effects jKnimholz & Dekelll20I ll : iKuhlen et alj|2012h . 
Moreover, if faint galaxies were the major contributors to reioniza- 
tion, this would imply an extended reionization epoch, a result dis- 
favoured by recen t constraints on kinetic Sunyaev-Ze l'dovich effect 
jZahn et al.|[2oTll : iKuhlen & Faucher-Giguerdl2012h . Therefore, if 
these dwarfs would not be able to provide the required contribu- 
tion to the ionizing background, additional assumptions have to be 
tested, as a strong luminosity/redshift evolution of /esc or the in- 
clusion additional sources of ionizing photons. 

All these considerations point out the need of a deeper investi- 
gation of the (possibly sub-dominant) contribution of AGNs to the 
ionizing flux, in order to explore to which extent these sources may 
be able to alleviate the tensions arising from considering galaxies as 
the only contributors to cosmic reionization. In this paper, we thus 
revisit and compare the relative contribution of AGNs and galaxies 
to the ionizing background at 2 > 5. In particular, we focus on 
the putative role played by faint AGNs, up to the same magnitude 
limits considered for galaxies. Both the determination of the faint 
end of the observed AGN luminosity function and the reliability of 
its extrapolation to fainter magnitudes present similar challenges 
with respect to analogous statistical estimators for galaxies. An ad- 
ditional source of uncertainty, however, arises from the apparently 
discrepant results obtained using independent selection techniques 
for the reference sample, i.e. X-ray or optically based, which show 
a different sensitivity to the various AGN populations. Nonetheless, 
it is also possible to combine the most recent determinatio n the 
galaxy stellar mass function of galaxies ( ISantinietani20I2l) with 
empirical arguments, to get an estimate of the maximum contribu- 
tion of the AGN population to the cosmic reionization. 

The structure of this paper is as follows: in sec. |2]we present 



the formalism we adopt to estimate the contribution of each class 
of sources with respect to the required ionization photon rate; while 
in sec. [3] we discuss the implication of our finding and in sec. |4] 
we present our conclusions. Throughout this paper we assume that 
QSOs represent the luminous sub-population of the homogeneous 
AGN population, and that Lyman Break Galaxies are a good tracer 
of the overall galactic population. 



2 MODELING THE CONTRIBUTION TO 

REIONIZATION OF DIFFERENT ASTROPHYSICAL 
SOURCES 

2.1 High-z QSO/AGN Luminosity Functions 

In order to estimate the contribution of quasars (QSOs) to the 
reionization we consider different observational constraints. A first 
set includes direct measurements of their high-z LF and its red- 
shift evolution. In this class, we consider both the optical QSO- 
LF derived in the framework of the Great Observatories Origins 
Deep Survey (GOODS) collaboration by iFontanot et"al] (l2007l 
F07 h ereafter) and the X-ray selected QSO-LF by Fiore et all 
(I2OI2I). Moreover, we also consider the upped limit estimate from 
IShankar & Mathud l l2007l) : they use the results of optical surveys 
down to the highest redshifts and faintest magnitudes proven (crit- 
ically including non-detections) to give constraints on the faint end 
of the QSO-LF at 2: > 5. Their analysis shows that the obser- 
vational constraints found so far are compatible with a faint-end 
slope of a = —2.8 (at 99% confidence level) and a = —2.2 
(at 90% confidence level). We directly compare these three esti- 
mates in Fig.[T] X-ray measurements and the relative analytical fits 
have been converted into absolute magnitu de A/uv by means o f 
the QSO bolomet r ic corrections propose d bv lMarconi et al.l(l2004h : 
lElvis et alJ ( ll994l) : lFontanot et alj ( l2007l) . The three estimates agree 
well at the bright-end of the LF; on the other hand the agreement is 
considerably reduced at the faint-end of the LF and at higher red- 
shift due to the different reds hift evolution extrapolated to higher 
redshift jFontanot et al.ll2007h . This behaviour is in part related to 
the different selection of QSO databases: it is indeed quite clear 
that X-ray surveys observe higher faint QSO space densities than 
optical surveys. It is also worth stressing that we expect obscura- 
tion effects to play a relevant role at these luminosities (see e.g. 
ISimpsonlliooi) . Moreover, at these luminosities we are no longer 
considering only "bona fide" QSOs (i.e. Mb < —23), but we are 
entering a regime where the AGN/QSO distinction becomes loose. 

If the QSO population and its high-z LF is representative only 
of a sub-population of the global AGN population, we may as well 
be missing a relevant contribution to the ionizing flux. As an al- 
ternative approach to gain insight into the AGN LF and in partic- 
ular its faint end slope we consider the galaxy stellar mass func- 
tio n (GSMF) as estimated in the redshift range 0.4 < 2 < 4.5 
bv lSantinietaLl ( l2012h . Following simple e mpirical reason i ng we 
tried to use the constraints provided by the ISantini et al.l ( 1201 2h 
analysis in their highest redshift bin (3.5 < 2 < 4.5). In partic- 
ular we assume: 



(a) each galaxy host a Supermassive Black Hole (SMBH) at its 
center, and the relation between galaxy s tellar mass and SMBH 
mass is the same as in local Universe ( lMagorrianetal.lll998l; 
iHiiring & Ri}3l2004 i.e. all galaxies are ellipticals); 

(b) the estimated faint end of the GSMF is a good representation 
of the dwarf population; 

(c) the high-z QSO-LF is a rescaled version of the GSMF; 
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Figure 1. Comparison of different estimates for tlie liigli-z QSO-LF. Blue 
lines and liglit blue shaded region correspond to the best fit models in 
iFontanot et alj )2007l) . and their redshift evolution. Green lines an d ligh t 
green shading coiTespond to the models in Shankar & Mathurl fcOOTt) . 
Yellow line corre sponds to the evolution of X-ray selected QSO-LF in 
iFiore et al ] i2012l) . Red lines and hatched light red shading correspond to 
the GSMF-derived high-z AGN-LF (see text for more details on the conver- 
sion fr om GSMF to Q SO-LF). Dots refe r to datapoints fromlFon tanot et al] 



bOOTl grey ast erisksXiFiore et al. I1 I2OI2I open circle s'). [ Civano et al jj2011l 

L et alj 



open sq uares), iBrusa et '^l " ]2010l open diamonds), iGUkmanl 
stars>. iJiang et alj j2009l open triangles). 
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(d) at z > 4.5 the GSMF evolves as the F07 QSO-LF: i.e. an 
exponential pure density evolution proportional to e^^'^^'^+~'; 

(e) at each redshif t only 0.45% of SMB H shine at their Eddington 
luminosity (see e.g. lHaiman et al.ll2004l) : 

(f) conversion from bolometric luminosity to Muv using a stan- 
dard approach (see e.g., F07): we consider a bolometric cor rection 
factor for the B-band luminosity of 10.4 l lElvis et al]|l994h and a 
B- to (71/-magnitude shift of -0.48 mag. 

We show the result of these conversions in Fig.[T| the red lines (and 
hatched light red shading in bet ween) correspond to the upper and 
lower best-fitflin lSantini et alj fcOl j) in the range 3.5 < z < 4.5. 
It is interesting to stress that the knee of the "mass-derived" lu- 
minosity function results in a good agreement with respect to F07 
best-fit models, and the overall normalization is consistent with the 
densities obtained in QSO-LF studies. This is an important san- 
ity check for the proposed approach, since some of the factors (a) 
to (f) are degenerate and uncertain as well. In particular, there are 
several evidences in favour of a redshift evolution of both the ratio 



^ In particular, the upper envelope coo'esponds to the best fit to a dou- 
ble power yielding a normalization iog('I>*) = —4.84, a characteristic 
mass log{Mi,) = 11.81, an high-mass-end slope /3 = —6.38 and a 
low-mass-end slope a = —2.27 (table 3 in ISantini et alj|2012l) ; while 
the lower envelope refers to the best-fit to a Schechter function yielding 
;oq(^„) = -4. 12, log{Mi,) = 11.30 and a = -1.80 (table 1 in 
ISantini et al.l2012l) . 



between SMBH mass and the spheroidal host component mass, and 
on the AGN fraction/duty cycle (see e.g. lShank^l2009l : lFiore et al.l 
2012 and discussion herein). The bright-end of the AGN-LF is 
ste eper than any QSO -L F, while its faint-end i s steeper than both 
the lFiore et"al] j20I2l) and lFontanot et al.l(l2007t) estiin ates. but still 



in the range allowed by the lShankar & Mathuij | |2007|) analysis. 

In order to compute the QSO contribution to reionization, we 
compute the rate of emitted ionizing photons per unit comoving 
volume Fagn as a function of reds hift, foll owing the same form al- 
ism as in lShankar & Mathud ( l2007l see also lMadau et al.ll999t) : 



rAGN(z)[s ^MpC 



o'u—. 



(1) 



p^(z)[ergs"^Hz"^Mpc"'^] = y ${L, z)L^{L)dL (2) 

In the above equations, Li, is in erg s^^ Hz^^, hp represents the 
Planck's constant, i^h = 3.2 x 10^^ Hz is the frequency at the 
Lyman Limit (i.e. 912 A); Vup = 12.8 x 10^^ Hz is the usual up- 
per limit to the integration, since photons more energetic are pref- 
erentially absorbed by He atoms. In practice, we assume that the 
absorbing cross section for neutral hydrogen a,^ is unity between 
i^H and i^up, and zero outside this range. This crude approxima- 
tion gives us a good grasp of the maximum contribution of QSO to 
the reionization background. We also assume that all the ionizing 
photons associated with the AGN spectra contribute to the ionizing 
background (i.e. fcsc ~ 1). In the calculations, we assume a QSO 
spectral continuum of the form — u~' and we assume a slope 
7 = -1.76 blueward of the Ly^, line ( iTelfer et al .1120021) . In the 
following, we will compare Fagn with the required total ionizing 
photon rate per unit comoving volume Fion, using the formalism 
prop osed in Madau et al. . (1999il . rescaled to WMAP7 cosmology 
as in ( |Pawliketalj2009l) : 



rion(2)[s"'Mpc-^] = 0.027 k ( ^ 



l + z 
7 



0.0465 



(3) 



where we convert the critical star formation rate into a photon rate, 
by assuming that k = 10^^'^ s~^ L yC photons per M© yr~^ are 
produced (see e.g. IShull et al.ll2012l see also eq. |6]below) and C 
refers to the cl umping factor of the intergalactic med i um. E arly 
work (see e.g. iMadau et al.l 1 1 9991 : IShankar & Mathud |2007|) as- 
sumed a high value C = 30 for the clumping factor, following 
the results of numerical simulations by Xlnedin & Ostriker ( 1 99"^, 
and concluding that the space density of ionizing photons deducted 
by observations was in most cases insufficient to reionize the Uni- 
ver se and/or kept it ionized. More recent theoretical estimates (see 



Bolton & Haehnellll2007l : IPawlik et alhoOSl : iHaardt & Madaiil 



[2013 revised the value of the clump ing factor towards lower val- 
ues. In particular. IPawlik et al. I( l200i estimates C = 6 as adequate 
for gas with densities of the order of the critical density for the on- 
set of star formation, while finding an even lower C = 3 value 
for gas with overdensities ~ 100. It is also worth stressing that C 
is expected to be a dec reasing function of redshift. For example 
iHaardt & Madaul j20I2h propose the following fitting formula: 



C{z) 



14-43 X 



(4) 



derived for gas with overdensit ies ~ 100. Eq. |4] is fully consis- 



tent with the lPawlik et al.l ( |2009|) estimate on the same overdensity 
scale. Lower values for the clumping factor reduce considerably 
the number of ionizing photons required to keep the Universe ion- 
ized at z > 6, mitigating the requests on the observed astrophysi- 
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cal sources. In the following, we will assume a redshift dependent 
clumping factor as in eq. |4] 



2.2 High-z LBGs Luminosity Functions 

In order to estimate the number of ionizing photons produced by 
the galaxy population, we consider the high-z luminosity func- 
tio n of Lyman B r eak G alaxies (LBGs hereafter) as estimated 
by iBouwens et alj ( 1201 ll) . Following these authors we describe 
the LBG-LF in the redshift range 3.5 ^2^8 as an evolving 
Schechter function (i. e. whose parameters a re evolving with cos- 
mic time, see table 1 in lBouwens et alj201lh . We then compute the 
luminosity density puv) using eq.|2]and we conve rt it to an estimate 
of the star formation rate density (psfr) using the lHaardt & Madaul 
( I2OI2I) conversion factor: 



PSFR(2:)[M0yr ^Mpc = 



puv(z)[ergs ^Hz ^Mpc ^] 
1.05 X 1028 



(5) 



Major uncertainties affect the conversion between UV luminosity 
density and the psfr. First of all, UV observatio ns are severely 
affected by dust attenuation. iBouwens et al.l ( |2007|) estimated that 
z > 4 LBGs suffer lower attenuation levels that lower-z counter- 
parts, and they suggest that the psfr obtained from eq.|5]may be 
underestimated by a factor ~ 1.5. Moreover, the conversion factor 
itself critically depends on the details of the Stellar Population Syn- 
thesis modeling; in particular the chosen Initial Stellar Mass Func- 
tion play a relevant role, by determining the relative abundance of 
massive stars, the main contributors to UV fluxes. The value we 
use in eq. |5] has been computed assuming a Salpeter IMF. If we 
consider, i.e. a Kroupa IMF, its value is reduced by a factor ^ 1.5. 
In the following, we take into account all these sources of uncer- 
tainties by assuming eq. |5] as a good representation of the mean 
conversion and by defining a "maximal" and a "minimal" model 
by increasing and decreasing the resulting psfr by a factor 1.5, re- 
spectively. We then estimate the rate of ionizing photon production 
using the convers ion factor proposed f or a lo w-metallicity gas by 
IShull et all 1 I2OII see also.Madau etlfll 19991) : 



rLBG(z)[s ^MpC ^] = K fasc pSFn{z) 



(6) 



where /esc represents the escape fractiorQ i.e. the fraction of pro- 
duced ionizing photons which are able to escape the local envi- 
ronment and ionize the intergalactic medium. This parameter has 
a key importance in order to evaluate the contribution of LBGs 
to the ionizing background, but its very poor const rained. At z ~ 
3 — 4, proposed values range from low (/osc < 5%, Vanzella et al.l 
20I0bl). to relativ ely high values (/esc ^ 20%. IShaDlev et al.l2006l : 



Iwata et al.ll2009l) . 



3 DISCUSSION 

3.1 The AGN contribution to the ionizing bacliground 

First of all we consider the AGN contribution to the ionizing back- 
ground. At variance with previous analyses we don't fix Lmin in 




Ullm 



Figure 2. Estimated AGN contribution to reionization at different cos- 
mic epochs. Lines and colours as in Fig. [T] Shaded areas represent al- 
lowed contributions after Helium reionization has been taken into account 
(see text for more details). Dashed vertical lines repres ent the actual limit 
(Mtjy ~ -20.6) of faint QSO searches as defined in lShankar & Mathud 
)2007h . 



eq. I2I but we study how the ratio Fagn /Lion evolves as a func- 
tion of redshift and limiting magnitude. We also compute the AGN 
contribution to the reionization by applying the same approach to 
our GSMF-derived AGN-LF (sec. 12. It . and include these results in 
Fig. [2] (red lines and red hatched area). Following the same (a) to 
(f) prescriptions described in sec. 12.11 we estimate that the puta- 
tive AGN luminosity associated to a Jeans Mass of pristine ga!0 is 
roughly A/uv ~ —10. The formation of the ancestors of SMBHs 
at very high redshifts ha s been discussed by a number of authors 
(see e.g. lPetri et al.ll2012l and references herein), showing that both 
heavy (> 10^ Mq) and light SMBHs seeds are physically plausible 
at the same scales accessible for gas cooling and star formation. We 
further assume that all ionizing photons associated with the AGN 
spectra are available for ionizing the IGM and we show our results 
in Fig.[2l 

We find that the QSO contribution at the current observational 
limit (Muv ~ —20.6), is not negligible at z ~ 6, but still in- 
sufficient to provide the required rate of ionizing photons. The 
difference with respect t o analogous works in the literature (e.g. 
IShankar & Mathudl2007h , is to be ascribed to the lower dumpiness 
adopted here. 

At higher redshift the contribution of AGNs to the ionizing 
background decreases rapidly, becoming of the order of a few per- 
cent at z ~ 8 and negligible thereafter. To obtain a relevant AGN 
contribution (i.e. > 10%) at these redshifts a steep faint end is re- 



2 Most observational works focus on the relative escape fraction, i.e. the 
fraction of escaping LyC photons, relative to the fraction of escaping non- 
ionizing ultraviolet photons. Since the relative /esc takes into account the 
dust attenuation, it is then possible to convert between the two determina- 
tions. For the purposes of this work we just deal with the absolute /esc- 



^ The limiting mass M j for the collapse of a cloud of cold gas into stars is 
usually detennined by comparing its thennal and gravitational energy; by 
assuming typical values for the temperature and density we obtain Mj ~ 
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Figure 3. Estimated LBG contribution to reionization at different cosmic 
epoclis. Yellow, blue and green lines and hatched regions refer to the inte- 
grated contribution of LBGs for /esc = 1%; 5%; 20% respectively. Solid 
and dashed line refer to the "maximal" and "minimal" models (see text 
for more details). The grey shaded area corresponds to the estimated AGN 
contribution from the iFiore et all | |2012|) high-z AGN-LF shown in Fig. |2] 
Dashed vertical lines repre sent the actual limit ( M^jy ~ —18) of faint 
LBG searches as defined in lBouwens et aljj201lh . 



quired for the high- z AGN LP (steeper than the GSMF faint-end 
slope determined bv lSantini et"ZI ( l2012h ') and an integration to very 
low luminosity limits, under the hypothesis that very faint AGNs 
are characterized by the same properties of their bright counter- 
parts. Another possibility is that the evolution of the low-luminosity 
AGN population becomes slower with respect to the 3 <«< 5 es- 
timates. 

We have also imposed that our analysis be consistent 
with present constrai nts on the Hell reionization (see e.g. 
iFurlanetto & Ohll2008h . i.e. areas of the parameter space predict- 
ing too many photons with energies larger than 4 Ryd and there- 
fore producing Hell reionization at redshifts jz 5 are forbidden. 
To this end we assume the same clumping factor as adopted for 
Hydrogen reionization, a ratio between the number densities of He- 
lium and Hydrogen atoms equal to Y/A{1 — Y) (where Y — 0.258 
represents the cosmic fraction of Helium by mass), a recombination 
rate for Helium 6 times faster than for Hydrogen and we integrate 
our standard QSO spectral shape between 4 and 16 Ryd. Shaded ar- 
eas in fig.[2]thus represents the allowed contribution to the ionizing 
background according to this analysis. 



3.2 The LBG contribution to the ionizing background 

In Fig.[3]we present the estimated contribution of LBGs to cosmic 

reionization for different fcsc values. Also for this class of sources 

we integrate the contribution of sources down to Muv ~ —10. 

We directly compare these results with the contribution from the 

I ' ' ri 

AGN-LF, and in particular with the results of Fiore et al. (201_2, 
grey shaded area). We confirm the result of Pawlik et al.i(.2009,) that 



observed LBGs may account for the total required ionizing photon 
budget at z ~ 6 if /osc 20%; at higher redshifts, we have to inte- 
grate the LBG-LF to increasingly fainter limiting magnitudes (up to 
A/uv ~ — 10 at z ~ 9), in order to produce enough ionizing pho- 
ton to fully account for cosmic reionization. However, if /osc ~ 5% 
for the whole LBG population, the LBG contribution is not enough 
to account for the whole required ionizing photon rate at 2: ;i, 7, 
even if we extrapolate the LBG-LF up to the fainter magnitudes, 
and it becomes roughly compatible with that of AGNs for z ^ 7. 
In order to achieve high-z reionization with LBGs only, a substan- 
tial contribution to the ionizing background of sources fainter than 
the actual observational limit is required and/or a very different (i.e. 
increasing) /esc with respect to their bright counterparts. 

In order to test this hypothesis we impose a simple luminosity- 
dependent /esc scenario to our maximal model by defining /osc ~ 
for Muv < —18.00 and a linearly increasing /esc with increasing 
magnitude: 



/esc = min[l,?7 x (Muv + 18.00)] 



(7) 



Despite the lack of constraints on the distribution of /esc in different 
galaxy population, this exercise allows us to provide a qualitative 
estimate for the typical magnitude of objects responsible for the 
bulk of reionization and their expected escape fractions, if /esc is 
indeed a decreasing function of luminosity. This simple toy model 
predicts that complete reionization at 2; = 8 by LBGs alone is 
achieved at typical magnitudes Muv ~ —13.5, —15, —15.5 for 
77 = 0.1,0.2,0.3 respectively. The predicted corresponding es- 
cape fractions are /esc 50, 60, 80%. A similar result is obtained 
at z — 9 (typical magnitudes Muv ~ -12.5,-14,-15 and 
/esc ''^ 55, 75, 90%). This results ar e compatible with the recent 
hydrodynamical simulations (see e.g.. lCiardi et al I20T I]), which re- 
quire rather high /osc values for fainter LBGs. It is also worth not- 
ing that, in general, lower /esc values are still compatible with an 
LBG-driven reionization if a redshift increasing efficiency of ion- 
izin g photons production and/or a top-heavy IMF are assumed (see 
e.g. Ischneideretal.ll2002h . 



3.3 Combined AGN-LBG contribution to ionizing 
background 

Fig.|4]shows the combined AGN-l-LBG contribution to the reioniza- 
tion of the Universe. For the sake of simplicity we only consider the 
maxi mum AGN contribu tion as estimated on the basis of the AGN- 
LF o f lFioreetal.l ( l2012l marked as a solid line in each panel): this 
is at the same time a conservative but representative choice among 
the various LF estimates that we consider in sec. 12. II From a com- 
parison with the results shown in Fig.|3]it is apparent that AGNs 
may significantly help reducing the gap between the LBG ionizing 
photon production rate and the required amount of ionizing pho- 
tons up to 2; ~ 7. For a given /osc the minimum luminosity of the 
galaxies required to match the theoretically estimated ionization 
limit is significantly increased. Nonetheless, at the highest redshift 
considered, it is still not possible to reach the expected space den- 
sity of ionizing photons for reionization if /esc 5%. Again, in 
order to reionize the Universe with LBGs and AGNs combined at 
such high redshifts, either /osc has to be higher than the cuiTent 
estimates (at least for the faint LBG population, which is assumed 
to provide the largest contribution), and/or the evolution of the LF 
has to slow down with respect to the present estimates. In Fig. |5] 
we then show the Fagn/Flbg ratio as a function of AIw and 
redshift. For the AGN population we consider as representative the 



© 0000 RAS, MNRAS 000, 000-000 



6 Fontanot et al. 




Ullm 



Figure 4. Combined LBG and AGN contributions to reionization. Colours 
and lines as in Fig. [5] In each panel we consider the max imal AGN con- 
tribution from the high-z AGN-LF from iFiore et al] ( l2012i marked by the 
solid line in each panel, see text for more details). 



20 -10 




Ullm 
™UV 

Figure 5. Logarithmic ratio between the AG N and LBG contri butions to 
the reionization. As a reference for AGNs, the lFiore et al.lll2012h AGN-LF 
has been considered. Colours and lines as in Fig. |3] The red filled triangle 
marks the position of our measurement at 3.4 < 2 < 4.0 canied out in the 
GOODS fields (see text for more details). 



-0.1 0.1 0.2 

Flux(AB) xl0-30[ergs-»cm-2Hz-'] 

Figure 6. Distribution of U-band fluxes (probing the rest-frame LyC) for 
67 LBGs (upper panel) and 7 AGNs (lower panel) in the GOODS South 
field (selected with 3.4 < 2 < 4 and 23.5 < 1775 < 26). Fluxes have 
been measured in both cases within an aperture of 1.2" diameter. The red 
dashed Gaussian distributions show the expected distribution of null detec- 
tions normahzed to the total number of measurements in each panel. For 
comparison, the dot-dashed blue histogram shows the expected distribution 
of fluxes fo r galaxies with an /esc = 10% (the IGM absorption being sim- 
ulated as in IVanzella et aDl2010bh . The vertical dashed line marks the 3 cr 
confidence limit. 



iFiore et al ]( l2012h LF, while for the LBG population we show re- 
sults for /esc = 1%, 5% and 20%. From Fig.|2]is clear that Tagn is 
decreasing at increasing redshifts at a faster pace relative to Flbg 
independently on f^sc- The latter contribution becomes dominant 
for faint sources at z > 5 for /esc ~ 5%. 

An interesting reference point can be obtained by studying the 
LBGs and AGNs in the GOODS South field, selected in equal red- 
shift and magnitude intervals, 3.4 < z < 4 and 23.5 < ^775 < 26. 
Following the procedure described in IVanzella et al. 1 ( l2010bh . we 
have measured the flux of 67 LBGs and 7 AGN in the U-band 
(probing the rest-frame LyC) using a circular apertur e of 1 . 2" diam- 
eter. A s shown in Fig.|6] one galaxy (dubbed lonl in lVanzella et al.l 
and two AGNs are detected above the 3 a confidence level. 
The average flux from AGN turns out to be (0.041 ± 0.08) x 10"^" 
erg~^ cm~^ Hz~^. In order to compute the corresponding quantity 
for galaxies we have considered that the distribution of their UV 
fluxes is characterized by 66 non-detections and one outlier. We 
have therefore assumed as an average UV flux of galaxies the flux 
of the outlier divided by 67 (we would obtain a similar value by 
averaging over the whole distribution) and as 1 a confidence lev- 
els those computed by Gehrels ( 1986) for small numbers of events 
(in this case one), obtaining (0.0024l^;[;°^|;) x 10"^° erg"^ cm"^ 
Hz"'^. The statistics are poor but allow us to roughly estimate a 
ratio between the AGN and LBG contribution to the UV ionizing 
background (shown in Fig. |5] with a red triangle) of 17^ {^3''. This 
value is conveniently independent on the IGM transmission and is 
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consisten t with the ^ 5% estim ate of the /esc from galaxies dis- 
cussed in IVanzellaeta"i] ( l2010bl) . 

It is worth noting that a simil ar measurement an d result has 
been obtained at lower redshift by ICowie et alj ( l2009h . These au- 
thors study the ionizing fluxes associated with the AGN and galaxy 
populations at z ~ 1.15 in the GOODS-North field by means of 
observations with the Galaxy Evolution Explorer (GALEX). Their 
results show the presence of a detectable signal corresponding to 
known AGNs/QSOs, while stacking analysis of galaxy images pro- 
vides no evidence for a significant ionizing flux (compatible with 



u. 



1%). 



4 CONCLUSIONS 

We have critically discussed in view of recent results the con- 
tribution to co smic reionization at z > 5 of both high-z QSOs 
jFontanot et al., 2 007: Shankar & Mathur 2007; Fiore et al. 2 oTl) 
and LBGs jSouwens et alJ l201ll) . In order to take into account 
the contribution from AGN fainter than the current observational 
depths we have also used a derivation of the A GN LF based on th e 
evolution of the galaxy stellar mass function (Santini et al.l l2012h . 
In the following we assume z ~ 7 (e.g. Mitra et al. 201^ as the 
fiducial redshift for a rapid transition of the hydrogen from a sig- 
nificantly neutral condition to a neutral fraction xhi « 10^'^. 
Our results show that 

(i) In order to achieve the HI reionization at z ~ 7 the properties 
of the AGN population have to be pushed to rather extreme values 
in terms of steepness of the faint end of the LF (a ^ — 2) and con- 
tribution of very faint (up to Muv ^ —10) objects. But in the case 
such conditions were met we would expect the reionization of Hell 
to take place, owing to the typical AGN SED, above r edshift z '--^ 5, 
which i s in con trast with present observations (e.g. iFechner et al.l 
( l2006h : lzheng et al. (2008)). AGN alone, at least in their standard 
manifestation, cannot be responsible for the reionization of HI. 

(ii) The LEG population may account for the whole photon 
budget needed for reionization, but only if the mean escape frac- 
tion of this population is of the order of /osc 20% and very 
faint dwarf galaxies provide a substantial contribution. Despite 
such an high fcsc is still compatible w ith observational constraints 
dShaplev et al.l2006l : llwata et al.l2009l) . several evidences point out 
that these values are not ty pical for the whole high-z LEG popula- 
tion l lVanzella et al .l2010bl seee.g.). If mean /osc valu es are indeed 
of the order of 5%, and more similar to z ~ 1 results jCowie et al.l 
l2009l) . the contribution of LEGs alone is not enough to account for 
cosmic reionization at z ^ 7, and we are forced either to advocate 
additional ionizing sources or a strong redshift/luminosity evolu- 
tion of /esc. 

(iii) If /esc ~ 5%, the AGN population can provide a significant 
contribution to the total photon budget and help achieving reioniza- 
tion not earlier than z ~ 7. 

(iv) If /osc ~ 5% for the brighter (L ^ L*) LEG galaxies a vi- 
able solution is to advocate an increasing /esc with decreasing lumi- 
nosity, as it may be expected if shallower potenti al wells are more 
easily cleared of (n eutral) hydrogen gas (see e.g. iMori et ^1 20021 ; 
lYaiimaetal .11201 ll) . A simple toy model with a linear growth of 
/esc with increasing magnitude shows that, in the case of a negligi- 
ble contribution to the ionizing rate from brighter LEGs, we need to 
assume /esc of the order of 70% for sources fainter than the present 
observational limits to provide a substantial contribution to reion- 
iz ation. These results are fully consi stent with the recent findings 
of lKuhlen & Faucher-Giguerg | |2012|) : these authors use a different 
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Figure 7. Comparison between observed and required ionizing background. 
Observed data from Bolton & Haehnelt (20071 circles), Becker et al. (2007, 
triangles), iFaucher-Giguere et al.l l i2008l squares), ICalverlev et al.l | |201 1 , 
pentagons). The black solid line mark the assumed minimum rate of emit- 
ted ionizing photons per unit comoving volume for Universe ionization. 
The cyan asterisk and dot-dashed lin e refers to the level c orresponding 
to the maximum contribution from the iFontanot et al.l l l2007h QSO-LF, in- 
tegrated up to Afiim = —10; the ye llow asterisk and dot-dashed line 
to the maximuin contribution from the iFiore et al.l ( 1201 2h QSO-LF, inte- 
grated up to Afiinj = —10; the dark green asterisk and dot-dashed line to 
maximum contribution from the Shankar & Mathur ( 20^^ QSO-LF, inte- 
grated up to Miini = —20; the dark green asterisk and dot-dashed line to 
maximum contribution from the GSMF-derived AGN-LF, integrated up to 
A7iim = —20; the blue asteris k and solid line to maximum contribution 
froin the iBouwens et al] )201 ih LBG-LF, integrated up to My^^-^ = — 10 
assuming /esc = 5%; the green asterisk and solid line to maximuin con- 
tribution from the LBG-LF, integrated up to = —18 assuming 
/osc = 20%; the red asterisk and solid line to the luminosity-dependent 
/osc scenario {rj = 0.1) integrated up to Miim = —14. The different slope 
for the two LBG-LF based models are due to the redshift d ependent evolu- 
tion of the LF shape, as inodeled bv lBouwens et all i201 ih results. 



approach to compute the contribution of LEGs to the ionizing pho- 
tons, based on the comparison of the LEG-LF with the expected 
Thompson optical depth and observational constraints on the Lya 
forest, favouring a scenario where /esc increases from z = 4 to 

z = 9. 

In Fig. |7] we compare our empirical predictions with the 
present observational constraints on the measured ionizing back- 
ground. We consider photoionization r ate data derived from the 
observed Lyg-forest effective opacity ( Eolton & Haehnell 1 20071 ; 
lEecker et al.l2007l;lFaucher-GiguWe et al.l2008h and from the QSO 
proximity effect ( Calverley et al.ll2oTl ) and convert them into ion- 
izing photons per unit comoving vol ume (Fbkg), following the 



same p rocedure and assumptions as in lKuhlen & Faucher-Giguera 
2|, their eq. 12), and using the me a n free path of ionizing pho- 



(1201 



tons as measured by jProchaska et al.l ( l2009h . The assumed mini- 
mum rate of emitted ionizing photons per unit comoving volume 
for Universe ionization below redshift 5 is only about a factor 2- 
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3 below the measured UV background. In this way, if the sources 
responsible for reionization at z ~ 7 have a rapid density evolu- 
tion, they would quickly exceed the constraint of the measured UV 
background. 

Fig. [7] show s that the space density evolution of AGNs 
(as estimated by iFontanot et"al] |2007| : IShankar & Mathurl |2007|: 
iFiore et al]|2012l) is too strong for these objects to provide a sig- 
nificant contribution to reionization 3.1 z ^7; LBGs show a milder 
density evolution for faint Mum values, steepening as /esc in- 
creases (and Afiim brightens), the steepeni ng beeing due to the re d- 
shift dependent evolution of the LF shape jSouwens et alj2()Ill) . 

In both cases our empirical models do not predict enough ion- 
izing photons at 2 > 7, if we force them to obey the constraints on 
Tbkg and /esc is kept constant. Also in this case, the most likely 
solution for an LBG-driven early reionization requires an increase 
of /esc either with decreasing luminosity (or with increasing red- 
shift): this is clearly shown by the red line in fig.|7] which repre- 
sents the contribution to the ionizing background in the luminosity- 
dependent /esc scenario (eq.|7] 77 = 0.1). 

The exact relative contributions from AGNs and LBGs criti- 
cally depends on the details of the slopes of the corresponding faint- 
end LFs, on the faintest luminosity limits of these populations and 
on /esc ■ Given the present uncertainties in its determination, better 
constraints on the high-z AGN-LF faint-end slope would be of fun- 
damental importance to understand the maximal contribution of the 
AGN population to the ionizing background. In general, our anal- 
ysis suggests that pushing the actual observational limits at least 
one magnitude fainter, despite very demanding from the observa- 
tional side, would be quite rewarding to clearly understand the rel- 
ative importance of different astrophysical sources in determining 
the ionization state of the early Universe. 
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